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We study theoretically the coupled multiferroic dynamics of one-dimensional ferroelec- 
tric/ferromagnet chains driven by harmonic magnetic and electric fields as a function of the chain 
length. A linear magnetoelectric coupling is dominated by the spin-polarized screening charge at the 
interface. We performed Monte-Carlo simulations and calculations based on the coupled Landau- 
Lifshitz-Gilbert and Landau-Khalatnikov equations showing that the net magnetization and the 
total polarization of thin heterostructures, i.e. with up to ten ferroelectric and ferromagnetic sites 
counted from the interface, can be completely reversed by external electric and magnetic fields, 
respectively. However, for larger system solely a certain magnetoelectrical control can be achieved. 
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I. INTRODUCTION 



Muttif envies, i.e. materials in which ferroelectric (FE) and ferromagnetic (FM) orders coexist, have attracted 
considerable attention in recent years [l|. The strong magnetoelectric (ME) coupling in these materials renders 
possible manipulation by external magnetic and electric fields 2-6]. This brings about entirely new concepts in the 
design of next generation devices. Currently studied multiferroics are of two kinds: single-phase systems and two- 
phase systems. Single phase multiferroics, such as RMn03 (e.g. R=Tb, Dy) [H, 0] have successfully been exploited 
to control the magnetic order by electrical means and vice versa. A number of interesting phenomena such as an 
electrically controllable persistent spin current Q and a modified Datta-Das spin-field-effect transistor |j| associated 
with the two dimensional electron gas formed at the interface of multiferroic oxides was reported. However, the small 
electric polarization and the low transition temperature in single-phase multiferroics [7j are adverse for applications. 
On the other hand, two-phase multiferroics [H, |T(|, including artificially synthesized ferroelectric and ferromagnetic 
materials, are very promising as a room-temperature ME system in which large ferroic orders coexist. 

Very recently, the dynamic response to applied external fields in a multiferroic chain with a linear ME coupling 
induced by electrostatic screening at the ferroelectic/ferromagnet interface was theoretically investigated It was 
demonstrated that for material parameters corresponding to BaTiOa/Fe the total polarization and the net magnetiza- 
tion are controllable by external magnetic and electric fields, respectively. In this work we study the size-dependence of 
the multiferroic dynamics of the electric polarization and of the magnetization in several heterostructures. Calculations 
are performed within the framework of the coupled Landau-Khalatnikov 12] and zero Kelvin Landau-Lifshitz-Gilbert 
equations [Hil ]. We find that the polarization (magnetization) are magnetically (electrically) switchable for chains 
with a small size. 



II. THEORETICAL FORMULATION 



As sketched in Fig[TJ we consider a two-phase multiferroic chain consisting of an artificially grown ferroelectric- 
ferromagnetic heterostructure such as BaTiOs/Fe. This system is similar to that in which the magnetoelectric effect 
was recently predicted theoretically [ill [ Ti l Il5j and realized experimentally [l6j . The total energy of the FE/FM-system 
in a general case falls into three parts [IJJ 
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The ferroelectric energy contribution reads 
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and the ferromagnetic energy part is 
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where E(i) and B(i) are the time-dependent applied external electric and magnetic fields, respectively. Various 
parameters and the values of reduced p^- and Si entering eqs. © and §5$ are discussed in details in Ref. in 
general there are several sources for the ME coupling. Here the ME coupling arises due to the accumulation of spin- 
polarized electrons at the FE/FM interface when the FE subsystem is polarized [17]. Since the electric field does not 
penetrate into the bulk of FM metals and the induced electric charges decay exponentially away from the interface, 
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FIG. 1. Multiferroic chain consisting of Nfe FE polarization sites and Nfm localized FM moments. 
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FIG. 2. The hysteresis loops of the reduced total polarization/magnetization of different ME chains in the presence of a 
harmonic external magnetic field B z (t) = BozCostot (left) or electric field E z (t) = Eozcosojt (right). The parameters are 
chosen such, that A = J, lj = 3.61 x 10 11 s _1 , olfm = 0.5. Amplitudes of external fields expressed in the energy units are 
Boz = J and Eq z = (left) and Bq z = and Eq z = J (right). 



the ME effect is confined to a depth on the order of atomic dimensions from the surface. Here, we consider an ideal 
FM metal, the exchange interaction between the induced surface magnetization and the FM part is thus limited to 
the first site only. The linear ME-coupling energy is given as 



E c = A Pl • Si 



(4) 



with A being the coupling strength. 

The magnetization dynamics proceeds according to the zero Kelvin Landau-Lifshitz-Gilbert[13| (LLG) equation of 
motion 
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where 7 = 1.76 • 10 11 (T s) 1 is a gyromagnetic ratio, «fm a Gilbert damping parameter. The total effective field 



1 5F 1 



acting on Si is defined as (t) ■■ 
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The polarization dynamics is governed by the Landau-Khalatnikov (LKh) equation as it was considered in the literature 
12] for ferroelectric materials 
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where 7^ is the viscosity constant. Numerically both eqs. ([5]) and ([B]) are solved in reduced units, renormalizing 
the energy (fTJ) over double anisotropy strength 2D. Then, the dimensionless time in both equations is r = uja t = 



7 B^t = 7 2D tj fig and the reduced effective fields are hf M (r) 
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III. RESULTS AND DISCUSSIONS 



The multiferroic response to external fields is investigated in two cases: (I) FE/FM chains with different lengths, 
Nfe — Nfm — 2, 20, and 200; (II) the first rii (n,- = 1, 2 and 20) polarization sites and magnetic moments around 
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FIG. 3. The hysteresis loops of the renormalized polarization/magnetization of the first n; (rn —1, 2 or 20) FE and FM sites 
as a function of a harmonic external magnetic field B z (t) (left) or electric field E z (t) (right) in the ME chain consisting of 200 
polarization sites and 200 magnetic moments. The parameters are chosen to be the same as in Fig. [TTJ 

the interface of the long chain consisting of Nfe = Nfm = 200 sites in total. Temperature is set to zero. This 
ensures the FM-subsystem remains away from its superparamagnetic limit and the FE-subsystem to be far from its 
phase transition temperature at nearly T = 10 K for the whole time scale of simulations. 

Aiming at fast switching devices we apply harmonic fast oscillating magnetic and electric fields, as in Ref. [ll[. We 
stress, however, that the switching by a magnetic field is not due to the resonance absorption of the magnetic field, as 
described in Refs. 18]. The strong magnetic field induces changes of the potential barrier compared with the height 



In the case (I) with different sizes of heterostructures, the system is described via the reduced net magnetization 
Ss(t) and reduced total polarization p s (t) 



Fig. HI! shows the ME response of the total reduced polarization/magnetization to harmonic magnetic field B z (t) = 
Bq z cos tot and electric field E z (t) = Eo z cosujt, respectively. Clearly, the FE polarization (FM magnetization) is 
completely switched by the external electric (magnetic) field. Additionally, the hysteresis loops indirectly driven by 
the applied field due to the ME coupling are presented. For the long ME chain, such as Nfe = Nfm = 200, the total 
FE polarization (magnetization) can not be reversed by the magnetic (electric) field since the active area of the ME 
effect is limited to the FE/FM interface only. A way around this problem is a multilayer stacking of alternating FE 
and FM structures coupled at their interfaces. 

In order to investigate whether only the several first polarization (spin) sites in the long ME chain follow the 
magnetic (electric) field, the first rii (n,= 1, 2 and 20) sites are extracted from the Nfe = Nfm — 200 chain. For 
convenience, we define the following renormalized parameters 
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The dynamical behavior is shown in Fig. 3. Although the total FE polarization (FM magnetization) driven indirectly 
by the external magnetic (electric) field, is not completely switched, it is indeed observed that the several first sites 
react strongly to the applied fields through the interface ME coupling. 
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IV. CONCLUSION 

In summary, the dynamical multiferroicity in a one-dimensional FE/FM heterostructure is investigated by means 
of the coupled LLG and LKh method as a function of the length of the heterostructure. These results are supported 
by our kinetic Monte-Carlo simulations, not shown here for brevity. We find that the FE polarization and the 
magnetization are switchable for thin MF-chains (e.g. Nfe = Nfm = 1,2 (Fig. |TTJ) and up to 10) by applied 
magnetic and electric fields, respectively. The complete MF-switching mediated by the multiferroic coupling is not 
achievable for larger MF-chains (e.g. for Nfe = Nfm = 20, as inferred from Fig. HU). Finally, we mention that there 
are several additional effects that may affect the multiferroic dynamics. In Ref. (TlJ , the effect of the depolarizing 
fields and the effect of the induced electric (magnetic) fields are discussed in some details. These effects are quite small 
for the BaTiOa/Fe chain and for the chosen range of frequencies. The temperature- and the frequency-dependence of 
the dynamical response deserves further considerations and is currently under study. 
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